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ABSTRACT 

We use fundamental plane (FP) distance estimates to the components of the double 
(S| ' cluster A2626 {cz ~ 17,500 km/s) to constrain cluster kinematics and estimate total 

O . binding mass. We employ deep R band CCD photometry, multi-object spectroscopy, 

^ ! and software designed to account for seeing effects to measure the FP parameters Re, 

O ! c, and (fie) for 24 known early type and SO cluster members. The FP coefficients from 

this sample {a = 1.30 ± 0.36 and f3 = 0.31 ± 0.06) are consistent with others reported 
in the literature. 

We examine the Mgh equivalent width distributions within both subclusters and 
find them to be indistinguishable. Lacking evidence for stellar population differences. 



■ we interpret the FP zeropoint offset between the two subclusters as a measure of the 

c/3 . distance difference. We find log^Ds/DA) = —0.037 ± 0.046, where Dd is the distance 

to subcluster d. This measurement is consistent with the subclusters being at the 
same distance, and it rules out the Hubble fiow hypothesis (distances proportional to 
^ ■ velocity) with 99% confidence; analysis of the subcluster galaxy magnitude distributions 

rules out Hubble fiow at 93% confidence. Both results favor a kinematic model where 
the subclusters are bound and infalling. 

We estimate the total cluster binding mass by modelling the subcluster merger as 
radial infall. The projected separation, the line of sight velocity difference and the line 
of sight separation constrain the cluster mass; the minimum possible total binding mass 
is 1.65 times higher than the sum of the standard virial masses, a difference statistically 
significant at the ~ 3a level. We discuss explanations for the inconsistency including 
(1) biases in the standard virial mass estimator, (2) biases in our radial infall mass 
estimate, and (3) mass beyond the virialized cluster region; if the standard virial mass 
is significantly in error, the cluster has an unusually high mass-to-light ratio (~ lOOO/i). 
Because observational signatures of departures from radial infall are absent, we explore 
the implications of mass beyond the virialized, core regions. 



^ Observations reported here were obtained at the MDM Observatory, a facihty jointly operated by the University of Michigan, Dartmouth 
College, and the Massachusetts Institute of Technology. 
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1. Introduction 

Galaxy cluster masses are clearly of cosmological 
significance, and examples abound: (f ) the differences 



in qluster evolution as a function of mass provide con- 



straints on the power spectrum ( Lacey fc Cole f993| . 



1994), (2) comparisons of the amounts of hmrinous 
and non-luminous matter in clusters provide clues 
to the nature of dark matter and the efhciency of 
galaxy formation ( David, Jones fc Forman 1995 ), and 
(3) the cluster mass-to-light ratio and baryon frac- 
tion constrain the cosmological density parameter Qq 
(e.g. iRamella, Geller & Huchra 1989|, [White et al 



1993, Mohr et al. 1996) 



Cluster mass studies tend to focus on the central 
regions where virial equilibrium is a more accurate 
approximation, and the weak lensing signals are eas- 
iest to detect. Here we apply a method to measure 
the mass beyond the virialized region in the double 
cluster A2626. We utilize recent improvements in the 
understanding of the fundamental plane (FP) to con- 
strain the total cluster binding mass by revisiting the 
subcluster radial infall model ( [Beers, Geller fc Huchra 
1982). Specifically, measuring FP distances (with 



an accuracy of 15% to 25%) to a reasonably large 
sample of cluster galaxies can yield a very accurate 



(sub)clustcr distance (Lucey & Carter 1988, Lynden 
Bell et al. 1988). For two bound subclusters, the 



line of sight separation, the line of sight velocity dif- 
ference, and the projected separation provide enough 
information to estimate the total cluster binding mass 
in cases where a radial infall model is appropriate. 

Although A2626 is rather distant, it is well suited 
for this analysis because simply demonstrating that 
the two subclusters are bound and infalling has inter- 
esting implications. A straightforward application of 
the radial infall model, using the (rest frame) line of 
sight velocity difference and the virial mass estimates 
within each subcluster, rules out bound and infalling 
merger models (Mohr, Geller fc Wegner 1996; here- 



after MGW96) ; the virial masses are not large enough 
to produce the apparent velocity difference. Thus, FP 
distances used to determine whether or not the sub- 
clusters in A2626 are bound and infalling serve as an 
independent test of the virial estimators, or alterna- 
tively, as a means of measuring the mass beyond the 
virialized region of the cluster. 

Possible pitfalls to this approach include (1) dis- 
tance biasing FP zeropoint (or shape) differences be- 
tween subclusters, (2) deviations from radial infall. 



and (3) significant interactions between the virial- 
ized regions of the two subclusters. FP variation 
with cluster environment is an area of active research 
(e.g. puzman et al. 1992 , Worthey, Trager, & Fabei 



T99|,]j0r gensen, Franx fc Kjaergaard 1996[; hereafter 



JFK96); current indications are that distance biases 
are smaller than 5% for ellipticals within clusters and 
are related to observable variations in the stellar pop- 
ulations within cluster galaxies. As discussed below, 
subcluster pairs suffering from pitfall numbers 2 and 3 
can be identified through corresponding observational 
signatures. We argue that the radial infall model and 
FP distances provide a promising method for deter- 
mining cluster masses on scales larger than the viri- 
alized region in low redshift clusters; if coupled with 
estimates of cluster light at larger radii, these mass 
measures should sufhce to determine whether mass 
to light ratios vary significantly between the virial- 
ized core and the surrounding infall region. For higher 
redshift cluster pairs, our approach could be used to 
quantify the level of environmentally induced funda- 
mental plane zeropoint variations. 

We describe the double cluster A2626, the FP ob- 
servations, and the reductions in §2. In Section 3, 
we derive FP coefficients, discuss constraints on the 
stellar populations in each subcluster, and derive con- 
straints on their relative distances. Section 4 con- 
tains a discussion of the radial infall model applied 
to A2626. Conclusions are summarized in Section 5. 
Throughout the paper we use Hq = lOO/i km/s/Mpc. 

2. Galaxy Sample and Data 



Below we describe Abell 2626 ( [Abell 1958[ ), sum- 
marize previous observational results detailed else- 
where (MGW96), present the new observations, and 
describe the reductions and analysis in detail. 

2.1. Abell 2626 



A redshift survey with the Decaspec (Fabricant & 



Hertz 1990[) mounted on the MDM Hihner 2.4 m tele- 



scope in FaU 1993 and 1994 revealed that Abell 2626 
((cz) 17, 500 km/s, richness class 0) is composed of 
at least two systems with mean line of sight velocities 
which differ by - 2, 500 km/s (MGW96). As is clear 
from Fig. |l] and Table |l|, there are two main compo- 
nents with cluster-like dispersions: group A centered 
at a velocity of (cz) = 16, 533 ± 141 km/s and group 
B centered at a velocity of (cz) = 19, 164 ± 138 km/s 
(90% statistical confidence limits). As discussed in 
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MGW96, the 11 galaxies centered at (cz) = 21, 173 ± 
119 km/s are most likely part of a background, low 
density structure because (1) the dispersion is low 
(cr = 200^52^ km/s) and (2) a large fraction of these 
galaxies are gas rich — 73% have emission lines com- 
pared to 39% (33%) for group A (B). If the velocity 
of the low density structure corresponds to Hubble 
flow then its distance from groups A and B is from 
20/1-1 Mpc to 50/i-iMpc. 

The group A and B masses, velocity dispersions, 
central densities and X-ray properties all differ. As 
noted in Table ^, the virial mass of subgroup A is 
roughly twice that of B, and the line of sight veloc- 
ity dispersions are a a — Q^S'^^]^ km/s and as = 
415^72^ km/s (90% statistical confidence limits). The 
absence of observed X~ray emission from group B can 
be used to place a limit on the central gas density 
in group B (MGW96). In particular, the Einstein 
imaging proportional counter (IPC) image of the re- 
gion reaches roughly 50 times fainter than the peak 
in the X-rav emission from group A. The group A 



cen tral gas densitv is ^ 3 x 10 cm ^. so the central 



gas i density in group B must be < 5 x 10 cm 



Alt hough the contrast in the projected galaxy densi- 



tie s of groups A and B is difficult to quantify, the 
smoothed distribution (Figure 3 in MGW96) indi- 
cates that the central density in group B is 7 times 
less than the central density in group A. Thus, the ra- 
tio of gas to galaxy density may be the same in both 
groups. 

As discussed below, if groups A and B are on a ra- 
dial infall trajectory, then the minimum implied grav- 
itational mass is more than the sum of their virial 
masses. However, it is also possible that the velocity 
difference is just due to Hubble flow. The magni- 
tude distributions of the confirmed members of both 
groups indicate that the merger hypothesis is favored 
over Hubble flow at 93% confidence (MGW96). To 
further investigate this issue we measure fundamen- 
tal plane distances to each cluster and determine the 
line of sight distance difference directly. To do so 
we first use R band CCD photometry to identify a 
sample of early type galaxies with redshift confirmed 
membership. 

2.2. Photometry 

On November 21, 1995, we used a thinned, Tek- 
tronix 10242 CCD mounted on the MDM Hihner 
2.4 m telescope to image galaxies in A2626 and J8 



( [Jackson 1982| ). The 24/z pixels are 0.275 " on a side. 
The seeing varied during the night between 0.85 " and 
1.05 ". Because the night was non-photometric, we 
determine the photometric zeropoint of each R band 
image externally. For the A2626 galaxies we use pho- 
tometric images taken previously (MGW96) on the 
MDM 1.3 m to zeropoint our images. Specifically, we 
use aperture photometry of isolated stars in the 2.4 m 
and 1.3 m frames to determine the 2.4 m zeropoint. 
The 1.3 m images are reduced to the Kron-Cousins 
system using Landolt (1992) standards. For J8 we 
rely primarily on p ublished photoelectr ic galaxy aper- 
ture photometry ( [Colless et al. 1993 ) to determine 
the R band zeropoint but also use stellar aperture 
photom etry from a photo metric 1.3 m image for one 
galaxy (Saglia et al. 1997). In addition to our primary 
imaging run, we acquired additional Hiltner 2.4m im- 
ages through the generosity of Paul Schechter; the 
reductions for these data were similar. 

The R band galactic extinction is 0.062 mag toward 
A2626 (NASA/I PAC Extragalactic Database, [Savage 



Mathis 1979) and 0.185 mag toward J8 (Saglia 



et al. 1997). We approximate the k correction as 



where Zg is the galaxy redshift ( Frei & Gunn 



1994). We account for the peculiar velocity compo- 



nent of the line of sight velocity when applying the 
cosmological dimming correction; specifically, two of 
the factors of (1 + z) in the cosmological dimming 
are due to relativistic effects and two are due to the 
change in the geometry of the universe between emis- 
sion and detection. Writing Zg = Zd + Zp, where z^ 
and Zp are the Hubble flow and peculiar velocity com- 
ponents of the redshift, the cosmological dimming is 
written C — -f 51og{l -I- Zd} + 5 log{l + Zg}; in terms 
of the observed and peculiar redshift, the correction 
is C = +10log{{l + Zg)y/1 - Zp/{1 + Zg)}. Thus, C 
is model dependent. As an example, if groups A and 
B are at the same distance and the peculiar velocity 



correction were incorrectly taken to be 



0, the 



fundamental plane distances would be biased so that 
group B would appear closer than group A by ~ 1%. 

We bias subtract, flat field (using twilight flats), 
and clean the images using standard IRAF tasks. For 
each galaxy, we calculate sky subtracted radial pro- 
files and uncertainties. Uncertainties have a Pois- 
son component added in quadrature with an assumed 
1% flat fielding uncertainty. We fit these profiles to 
summed -R^^'^laws and exponential disks. 



Ib{R) = hoe' 



-7.67[(i?/fl6)i/''-l] 



(2-1) 
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where I(,o (Rt) a-nd Ido (Rd) are the characteristic 
surface brightnesses (scale lengths). Our fitting pro- 
cedure accounts for the effects of PSF smoothing and 
pixel extent as described in [Saglia et al. 1993b . 



More specifically, we determine the parameters of 
the i?^''^and disk components by minimizing the 
difference between the seeing convolved model and 
the observed profile. The expected surface brightness 
I{R, dR) in the bin of width 2dR and radius R is 

rr/? ,j,. [F{R + dR)-F{R-dR)] 

where F{R) is the seeing convolved integral flux 
within radius R. The integral flux follows from the 
bulge, disk and PSF parameters. Specifically, 

F{R) = rJ^ dkJi{kR)p{k) (/b(fc) + /rf(fc)) (2-3) 

where Ji is a Bessel function, and p{k), Ib{k), and 
Id{k) are the Fourier transforms of the PSF, the bulge, 
and the exponential disk. Following Saglia et al 
19S(3^, we use a PSF of the form 



m 



(2-4) 



and find it to be a good fit to stars on the 2.4 m images 
with 1.4 < 7 < 1.6. We use the analytic approxima- 
tion to the Fourier transform of the i?^/^law ( ^aglia 
et ;.l. 1993b), and the exact form of the transform of 



the disk component (e.g. Bracewell 1986 ) 
The fitting procedure has three steps: 



• Fit the PSF parameters b and 7 using radial 
profiles of 4 stars within each 2.4 m image; 
best fit parameters minimize the difference 
between the observed and model profile. 

• Feed initial guesses for the bulge and disk pa- 
rameters into a s implex minimization routine 

which minimizes the dif- 



( Press et al. 1 



ference between the seeing convolved theoretical 
and observed profiles. At each iteration the full 
Fourier integral is solved for each point in the 
radial profile; with an HP-735 the minimization 
takes several minutes. In addition to the un- 
certainties in the observed profile, a theoretical 
uncertainty is introduced to account for imper- 
fections in the smoothing corrections (e.g. the 



exact form of the PSF) . Specifically, the uncer- 
tainty for each profile point has a component 
which is proportional to the fractional effect of 
the smoothing on the model profile. This has 
the positive effect of making the fit results less 
sensitive to the inclusion/exclusion of the heav- 
ily smoothed inner points in the profile and to 
the details of the smoothing operation. We in- 
clude a theoretical uncertainty which is 10% of 
the fractional effect of the smoothing on the 
model profile. This is consistent with a phi- 
losophy that we can correct for seeing effects at 
the 90% level. 

• We fit a bulge only model to the profile and 
compare the minimum x^ for this model to the 
bulge plus disk x^ . For galaxies where x^ is only 
minimally improved by the inclusion of the disk 
component we use the bulge only fit. Once the 
fitting is complete, we use the best fit parame- 
ters to calculate R^- the half luminosity radius, 
fif,— the surface brightness at i?e, and {fJ.e)~ the 
mean surface brightness within R,,. 

We test our fitting routines by (1) fitting artificial 
galaxy images, (2) comparing results from multiple 
images of the same galaxy, and (3) comparing our re- 
sults to those of Saglia et al. (1997) in the cluster 
J8. We create artificial images by adding a galaxy 
described by an R^^^ profile to a flat background of 
10'^ cts/pixel. We then smooth the artificial image 
and introduce Poisson noise. By using delta func- 
tions instead of i?^/^ profiles we produce stars which 
are then used to measure the PSF parameters b and 
7. The process of fitting artificial galaxy images tests 
the accuracy of our fitting in that case where the 
true galaxy profile is well described by an R^^^ pro- 
file. We find that the parameters Re and (/ig) are 
constrained to an accuracy of Alog(i?e) ^ 0.02 and 
A (fie) 0.07 for Re as smaU as 50% the FWHM 
of the PSF. As expected, the accuracy of the com- 
bination log(i?e) — 0.33 * (/ie) (which appears in the 
equation describing the fundamental plane) is con- 
strained to ~ 0.01 because the errors in Re and {fie) 
are correlated. 

We also compare results from subarcsecond seeing 
images of seven A2626 galaxies to results from poorer 
quality images. The mean offset between the param- 
eters in the good and poor seeing images of the same 
galaxies are consistent with zero, and the scatter is 
smafl. Specifically, A{fie) = -0.0262 ± 0.0311 with 
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an RMS of 0.0823, Alog(i?e) = -0.0073±0.0074 with 
an RMS of 0.0196, and A(log(i?e) - 0.33* (^fe)) = 
0.0013 ± 0.0037 with an RMS of 0.0097. 

We observed five galaxies in the EFAR cluster J8 
(Wegner et al. 1996) for which we have photometric 



zeropoints. These observations provide the strongest 
test of our methods because the data, reductions and 
analyses differ significantly. There are no significant 
offsets between our values and EFAR values (Saglia et 
al. |1997| ); specifically, A (/ie) = -0.0766±0. 1223 with 

-0.0257±0.0326with 

0.33*{^^e)) = 
These three 



an RMS of 0.2735, Alog(i?e) 
an RMS of 0.0729, and A (log(i?e 
-0.0006±0.0112with an RMS of 0.0249. 
tests of our photometric methods indicate that our 
procedure is accurate and makes only a small con- 
tribution to the observed scatter in the fundamental 
plane (see below). 

2.3. Spectroscopy 

On November 22 and 23, 1995, we used the De- 



caspec (Fabricant & Hertz 1990) plus Mklll spectro- 



graph mounted on the MDM Hiltner 2.4 m telescope 
to obtain spectra of four galaxy fields in Abell 2626, 
one field in J 8 ( Wegner et al. 1996| ), and multiple 



spectra of five HD stellar velocity standards. The De- 
caspec fibers subtend 2.3" in the focal plane. We used 
a g-psm blazed at 5.700 A with 600 1/mm, yielding a 



typlical PSF of 5.3 A at 2.18 A per pixel and coverage 



from 4,300 A to 6,500A. Four or five thirty minute 
exposures were obtained for each of the A2626 fields; 
shorter exposures were obtained for the brighter J8 
field. We extracted the spectra, using arc lamp ex- 
posures taken on either side of each object exposure 
to dispersion correct. We then combined the spec- 
tra from the four to five exposures, excluding cosmic 
rays through sigma clipping of large, positive devia- 
tions. With the Decaspec there are four sky spectra 
for each object spectrum. We combined the sky spec- 
tra for each object, scaled the sky spectrum using the 
flux in the 5577 A sky line, and then removed the 
sky contribution from the object spectrum. The final 
spectra range in signal to noise (S/N) from S/N=20 
to 75 per pixel at 5,300 A. 

We obtained spectra of more galaxie s, stellar ve- 
locity an d Lick linestrength standards ( Worthey et 
il. 1994) using the Modspec mounted on the MDM 
Hiltner 2.4m telescope in Nov and Dec '96. We used 
a 1,200 1/mm grating, the 200 mm camera, a thinned 
CCD with 24 fim pixels, and a 1.7" wide slit; stan- 
dard longslit reductions led to spectra with 2 A reso- 



lution, 1 A/pixel sampling, and coverage from 4,800A 
to 5,800A. Galaxy exposure times ranged from 20 min 
to 2.5 hrs, and the spectra have S/N> 30 per pixel. 
The higher resolution of this setup allowed us to accu- 
rately measure galaxy dispersions as low as 100 km/s. 

For each run, galaxy dispersions and velocities are 
extracted using the stellar templates from that run 
and the cross correlation peak fitting available in fx- 
cor (IRAF.RV); we use the approximate wavelength 
region 4,800-5,800A which contains the strong Mg 
features. The correlation peak width is transformed 
to an intrinsic galaxy dispersion using translation ta- 
bles created for each stellar template; the transla- 
tion tables are produced by cross correlating the stel- 
lar templates against Gaussian broadened versions of 
themselves. The final galaxy dispersion is the mean 
of the dispersions from each stellar template in the 
run. Finally, an aperture correction 



log 



fobs 



0.0381og(^^) (2-5) 



is applied. This aperture correction accounts for the 
falling dispersion profile in early type galaxies; we 



normalize to an aperture of radius 



ap 



1.7" and a 



galaxy with 9" = 20". For the longslit spectra we use 
fap = 1 .025 xy / TT , where x and y are the width and 



length of the rectangular aperture (J0rgensen, Franx 



fc Kjaergaard 1995| , [Baggley 1996|) 



We estimate the uncertainties in our dispersions 
using Monte Carlo techniques. For each galaxy spec- 
trum we broaden a stellar template to mimic the mea- 
sured dispersion, scale the template to the observed 
cts/pix in the region around 5,200 A, and add the sky 
at the observed level. We sample this artificial galaxy 
spectrum 100 times, introducing the Poisson noise 
and sky subtracting, and then measure the dispersion. 
We use the RMS of the measured dispersions around 
the input value as a measure of the velocity dispersion 
uncertainty (see Table H). Not surprisingly, the sim- 
ulations indicate that averaging the dispersions mea- 
sured by cross correlating against multiple templates 
does not improve the accuracy. However, these sim- 
ulations consider only Poisson noise effects; we aver- 
age the multiple dispersion measurements from each 
galaxy spectrum to reduce template mismatch sys- 
tematics. 

We have multiple spectra for 21 galaxies. The dis- 
tribution of dispersion differences scaled by the un- 
certainties for the 28 comparison pairs has an RMS 
of 0.97, indicating that the Monte Carlo uncertain- 
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ties are a reasonably good estimate of our true dis- 
persion uncertainties. The mean scaled difference 
{{{o'Mod — o'Deca)/£a-) wherc Ca is the Uncertainty in 
the difference) between Modspec and Decaspec mea- 
surements for 23 comparison pairs has a mean of 
—0.11 ± 0.19 with an RMS scatter about this mean 
of 0.89; the variance weighted, mean logarithmic dif- 
ference between Modspec and Decaspec dispersions 
is —0.001 ± 0.013. Because there is no evidence for 
a systematic difference between dispersions measured 
wit h the Decaspec and the Modspec, we apply no cor- 



sion correction we determine by broadening our Lick 
standards, and Cap is an aperture correction of the 
same form as in Eq. 



2-5 but with a coefficient of 



rcc pion. We note that the Modspec was used in mca- 
suring dispersions for 10 of the 16 galaxies in group A 
and 5 of the 8 galaxies in group B; so any systematic 
offset between Decaspec and Modspec measurements 
would add to the FP scatter within each group, but 
would not bias the estimate of the distance difference. 

Table || contains a list of the galaxies and their 
properties; for the galaxies with multiple measure- 
ments, we use the variance weighted average disper- 
sion. The columns of Table || contain the galaxy tag 
(first letter corresponds to group membership), the 
coordinates, the redshift cz, the velocity dispersion a 
and uncertainty Eg-, the R band apparent magnitude 
Mji, the mean surface brightness (/ig) within the half 
light radius, the half light radius 9^, the Mgi, equiva- 
lent width, and the bulge luminosity fraction Fb- 

For 6 galaxies in J8 we compare our measured dis- 
persions with those of Saglia et al. (1997). There is 
evidence of an offset; the distribution of {tefar — 
o'obs)/^a has a mean of — 0.80± 0.31 (where is the 
uncertainty in the dispersion difference). The aver- 
age EFAR velocity dispersion is 15 ± 7 km/s lower 
than ours. This offset does not affect our estimates 
of the relative distances to the two subclusters, but 
would complicate efforts to bring our distances onto 
the EFAR system. 

We place constraints on possible differences in the 
stellar populations of the two subclusters by using 
IRAF scripts to measure Mgt, linestrengths in each 
galaxy spectrum. The rest wavelengths given in 
Burstein et al. (1984) for the feature and continuum 
bands are corrected using the redshifts in Table |^. We 
broaden our spectra to the nominal 8.6 A resolu tion 
of the Lick system ( Worthey fc Ottaviani 1997|) be- 
fore making measurements. We tranform to the Lick 
system using an expression of the form 



0.050 instead of 0.038. We determine a separately 
for each run; observations of Lick standards indicate 
a = 1.060±0.016 for the Nov '96 and a = 1.086±0.028 
for the Dec '96 data. We use overlapping galaxy 
observations to measure a = 1.066 ± 0.034 for the 
Nov '95 data. The aperture correction is similar to 
that used by J0rgensen (1997), but generaliz ed to the 
form used for the velocity dispersions above ( Baggley 
1996|) . 



We test the accuracy of our Mg^, measurements us- 
ing multiple galaxy observations. Twenty three com- 
parison sets from 18 different galaxies yield an RMS 
of 0.051 in AlogMt;;,. Thus, we estimate the uncer- 
tainty of a single observation is approximately 9%. 
The linestrengths listed in Table ^ are averages of 
multiple observations where appropriate. Finally, we 
compare our linestrengths to preliminary EFAR Mgi, 
linestrengths in 6 J8 cluster galaxies; the EFAR values 
are lower: AlogM.gb = -0.0733 ±0.014 with an RMS 
of 0.039. The hnestrengths in Table | are corrected 
to the EFAR system. 

3. Fundamental Plane Analysis 

In this section we, determine the coefficients and 
zeropoint of the FP within the two clusters (§3.1), 
discuss the evidence for stellar population variations 
(§3.2), and then interpret the observed FP zeropoint 
offset (§3.3). 

3.1. Determining FP Coefficients 

Even if the assumption that elliptical and SO galax- 
ies within both clusters are similarly distributed within 
the FP is valid, the FP coefficients and cluster dis- 
tances can be biased by differences in selection (e.g. 



(2-6) 



where a is a scale factor, /(ct) is a velocity disper- 



Lynden-BeU et al. 1988| , [Baggley 1996|) . Our galaxy 
sample is drawn from a list of known elliptical and 
SO cluster members sorted by central aperture mag- 
nitude (aperture is 2" x 2" square) ; the parameters of 
our final sample indicate that central aperture mag- 
nitude is correlated with log^e — 0.3 (/^e), a quantity 
similar to the combination which appears in the FP. 
Thus, our selection is similar to selection by isopho- 
tal diameter Z?„ (Dressier et al. 1987). We seek to 
minimally bias the coefficients by (1) combining the 
samples from both clusters into a single FP fit and 
(2) weighting galaxies within each cluster by the sam- 
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pic completeness. Finally, we examine the variation 



of Ihe zcropoint offset between the two clusters as a 



function of variations in the FP coefficients (§3.3). 

Our galaxy sample extends over 1.5 mag in cen- 
tral aperture magnitude. We divide this range into 
four bins and calculate the completeness within each 
separately for the two clusters. The completeness in 
each bin is simply the number of galaxies within our 
final FP sample divided by the estimated total popu- 
lation of elliptical and SO cluster members. We esti- 
mate the total cluster population within each bin to 
be Nci = Ntot * fch where Ntot is the total number 
of elliptical and SO galaxies identified in the R band 
survey of the cluster, and f^i is the fraction of all el- 
liptical and SO galaxies with redshifts that lie within 
the velocity range of cluster cl. 

Weighting by the inverse of the completeness ap- 
pears to be ill advised for the faintest galaxies within 



orthogona l deviations from the plane (JFK96, Bag- 
gley 1996| ). The best fit values are a — 1.30 and 



(3 — 0.31; the scatter in logi?e is 0.090 (23% distance 
uncertainty). Our FP has the form 

log6ie = (1.30 ± 0.36) logcr + (0.31 ± 0.06) (/le) + 7- 

(3-2) 

The uncertainties are determin ed using a bo otstrap 
resampling procedure (JFK96, Baggley 1996 ) to re- 



sample the galaxy list and fit the FP 10 times; the 
uncertainty intervals are half the width of the 68% 
confidence region. For comparison, the FP coeffi- 
cients in the case where all galaxies are given equal 
weight are a = 1.27 and P = 0.27. Our best fit FP 
coefficients are statistically consistent with the val- 
ues a = 1.24 ± 0.07 and /3 = 0.33 ± 0.01 determined 
in a sample of 207 galaxies with Gunn r photome- 
try (JFK96); as discussed by JFK96, their coefficients 



eacii clusLtii. Fui Lhesu Lwu galaxies Llit; cunipluLenuss 



is s nail enough (wLiighL is 13.2 (17) fui Lhe faiiiLesL 



are consistent with previous studies (e.g. Faber et al 



Trr 



1987, Djorgovski & Davis 1987, Bender, Burstein & 



gro 



up A (D)), LhaL weigliLiug by Lhe iiicumpleleuess 



Faber 1992, Guzman, Lucey & Bower 1993, saglia. 



would lead Lu Lhem duiuiiiaLiug Lhe flL. We exclude 



the'faiiiLesL galaxy fiuiu each sample when deLeiiuiii- 
ing the FP coefficients, but include them when cal- 
culating the best estimate of the cluster zeropoint. 
Finally, we normalize the weights within each cluster 
so that the ratio of weighting for cluster A:B is 15:7, 
the ratio of the number of galaxies in each. 

We fit the galactic properties to a FP of the form 
log — alogtr -|- (3 {^e) + 7, where a and /? are 
the (cluster invariant) coefficients of the plane and 
7 is a distance dependent zeropoint (Faber et al. 
19 ^, [Djorgovski fc Davis 1987| , p0r gensen. Fraud 
& ijaergaard 1993, JFK96). We combine galaxies 
from both groups in a single fit by using deviations 
around median values of distance dependent quanti- 
ties (Baggley 1996). The corrected quantities {^e) 



Bender & Dressier 1993, J0rgensen, Franx fc Kjaer- 
gaard 1993)! Fig. ^ contains plots of the group A 



and logcr are distance independent, so we transform 
log^e alone; specifically, we remove an estimate of the 
zeropoint within each subcluster by expressing log Of, 
as the variation around the which minimizes the 
cluster scatter in 7 



(filled) and B (hollow) galaxies along two projections 
of the FP; the lines in each figure represent the FP 
for groups A (solid) and B (dashed). 

3.2. Constraints on the Stellar Populations 

We compare the distribution of Mgjj equivalent 
widths within the two subclusters to place constraints 
on stellar population variations which could introduce 
offsets in the cluster FP zeropoints, biasing the rela- 
tive distances to the two clusters (e.g. JFK96). Fig. || 
contains a plot of log {Mgt) versus log a for the galax- 
ies in groups A (solid points) and B (hollow points). 
The lines represent best fit models for groups A (solid) 
and B (dashed). 

Rather than allowing both the slope and zeropoint 
to vary, we constrain the slope to be that found by 
J0rgensen (1997) for a sample of ^ 300 galaxies. 
Specifically, the combined correlations of Mg2 with 
Mgh and Mg2 with a imply 



log 8e -> log 9e - 7/»t 



(3-1) 



logM^b = 0.3071ogcr-5, 



(3-3) 



Note that the transformation depends on the FP co- 
efficients, so it must be reapplied during each fit iter- 
ation. 

We determine the best fit FP coefficients by mini- 
mizing the sum of the absolute value of the weighted. 



where the zeropoint for the large galaxy sample is 
b = 0.034 (J0rgensen 1997). By minimizing orthogo- 
nal deviations from the line, we find b = 0.054± 0.016 
for sample A and b = 0.042 ± 0.022 for sample B. 
The zeropoint of each cluster is statistically consistent 
with the zeropoint from the J0rgensen sample. We 
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can measure the offset between the two clusters some- 
what more accurately; the zeropoint difference be- 
tween the two clusters is Ab = bA—hs = 0.012±0.019, 
statistically consistent with no offset. 

Trends in AIg2 have been noted as a function of 
radius from the center of the Coma cluster (Guzman 
et i J. 1992D and as a function of cluster velocity dis- 
persion for a sample of 11 clusters (JFK96; J0rgensen 
19S|^. Another study indicates that elliptical galax- 
ies in low velocity dispersion environments have both 
higher metallicities and younger stellar populations 
than galaxies in higher velocity dispersion environ- 
ments (Rose et al. 1994). The velocity disper- 
sion difference between clusters A (ct = 658 km/s) 
and B (cr = 415 km/s) should introduce an offset 
of Ab ^ —0.016 into the Mgi,~a correlation. This 
expectation is inconsistent with our observations at 
the 1.5(7 level; we would require additional galaxy 
linestrengths to resolve an offset of that magnitude. 

3.3. FP Zeropoint Differences 

We determine the zeropoint of each group by using 
the median 7; the 16 galaxies in group A yield ja = 
-8.4083 with an RMS around this value of 0.0879, 
and the 8 galaxies in group B yield 73 — —8.3632 
with an RMS of 0.1149. The homogeneous Malmquist 
bias has a marginal effect on the relative distances to 
these two clusters (Lynden-Bell et al. 1988); we cal- 
culate 6jA = —0.0033 and 5-fB — —0.0114, resulting 
in an increase in A7 of ~ 0.0081. Thus, the zeropoint 
difference between the two clusters is 



is the distance of galaxy x and F follows from the 
standard form implied by the FP 



A7 = 7A - 7B 



-0.037 ±0.046, 



(3-4) 



where we have assumed the uncertainty in the zero- 
point is the scatter about that value divided by \/N, 
where N is the number of galaxies (note that this is 
a conservative estimate of the true uncertainties be- 
cause the scatter is sensitive to outliers). As discussed 
in §2.2, if we assume groups A and B are at the same 
distance rather than at distances proportional to their 
mean redshifts, the cosmological dimming corrections 
change so that A7 increases by ~ 0.004. 

In general, zeropoint offsets can be caused by stel- 
lar population or metallicity induced M/Lji differ- 
ences and/or distance differences. Because there is no 
compelling evidence for stellar population differences, 
we use the observed zeropoint offset and uncertainty 
as relative distance constraints. The quantity 7 (in 
Eqn. |3-2| ) can be written — ^og{T/ Dx) where 



Re = Fa" {le 



-2.5/3 



(3-5) 



Assuming the intrinsic FP zeropoint log F is the same 
for all galaxies, the relative distances to two galaxies 
X and y is 

Fig. ^ contains a plot of the distances to the galaxies 
in groups A (solid) and B (hollow) relative to the me- 
dian distance to group A versus the galaxy redshift. 
Also marked with vertical lines are the (homogeneous 
Malmquist bias corrected) median distances to groups 
A {Da; solid line) and B {log {Db/ Da) = -0.037; 
dashed line); the horizontal lines mark the mean ve- 
locities for groups A (solid) and B (dashed). The error 
bars on the dashed, vertical line represents the rela- 
tive distance uncertainty. The large star marks the 
position of cluster B if its distance relative to cluster A 
were reflective of pure Hubble flow (AD (x A{v)); the 
Hubble flow hypothesis implies log {Db/ Da) = 0.065. 

The zeropoint difference is consistent with zero, 
and it is sufRcient to rule out Hubble flow at 2.2f7; 
therefore, the FP measurements suggest that cluster 
B is closer to us than cluster A (in §5.1 we return to 
this issue in more detail). Naturally, the zeropoint 
offset depends on the coefficients of the FP, which are 
only constrained to within ^20% of the best fit values. 
Using the distribution of 7^ — 7s from the 10^ boot- 
strap resampling fits, we find that in 99.05% of the 
simulations the zeropoint offset between the two clus- 
ters is less than the Hubble flow offset, rejecting the 
Hubble flow hypothesis. This probability corresponds 
to a -1-2.35(7 deviation in a Gaussian, similar to the 
2.2(7 estimate above which we derived from the scat- 
ter about the median distance in each cluster. The 
bootstrap probability is superior because it includes 
variations in the FP coefficients (and is less sensitive 
to outliers than the scatter estimator). 

4. Radial Infall Model as a Mass Constraint 

The radial infall model has been applied to many 
close galaxy cluster pairs to determine whether or not 



the pairs will merge (e.g. 


Beers, Geller & Huchra] 


1982, 


Beers et al. 1991, Beers et al. 1992, |Colless & 


Dunn 1995, 


Scodeggio et al. 1995|). In the following 



we use the measured line of sight separation and this 
model to estimate the total binding mass required to 
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explain the observed infall velocity. We then compare 
this mass to virial estimates of the cluster binding 
mass. 

4.1. Radial Infall Model 

For two isolated, bound objects within an expand- 
ing universe, the dynamics can be parametrized in the 
standard way 



t =^ B(ri — sin 77) 
1= (GMoB'^y^^ {1 - COST]) 



(4-1) 



V B 



sm?7 



(1 — cos 77) 

where ry is a development angle, I is the separation, 
V is the relative velocity, t is the age of the universe, 
B is an undetermined coefhcient, and Mg is the sum 
of the individual masses (e.g. Peebles 1993). When 



we observe two clusters on a merger trajectory at the 



we term the cluster merger mass Mq). Mq can then 
be compared to other mass estimates, which typically 
probe the cores of clusters. 

Departures from radial infall and the physical ex- 
tent of the clusters can bias the merger mass. Gen- 
erally speaking, the large mass and rarity of galaxy 
clusters make them good candidates for this analysis. 
We plan to study the accuracy of merger masses us- 
ing numerical simulations of cluster evolution within 
"realistic" environments. 

Note that there are scenarios in which the merger 
mass would be substantially biased; it is possible to 
recognize these scenarios observationally. In the first 
scenario, the two merging clusters are at small sep- 
aration (< 1 Mpc); interactions among the cluster 
components, significant overlap of their mass profiles, 
and small departures from pure radial infall combine 
to bias the inferred merger velocity and the implied 
merger mass. Fortunately, the signatures of cluster 



present epocti, we measure the projected separation 
of ihe clusters Al± and a ime ot sight velocity dit- 



mergers have been extensively studied (e.g. Geller & 
Beers 1982 , Dressier fc Shectman 198§|, Jones fc For- 



fersnce Avj^qs- Modelling the collision as a radial 



man 1992, Mohr, Fabricant fc Geller 1993, Pearce 



Thomas & Couchman 1994, Mohr et al. 1995 



infill, we parametrize the merger m terms ot the an- 



gle between the merger trajectory and the line ot 
sight {Avlos = Aw cos and AZ_i_ = A/sin0). 

With four free parameters in this model (jj, B, cf) 
and Mg) and three observables (the age of the uni- 
verse tn, the projected separation A^j^ and the line 
of sight velocity difference Avlos), the total cluster 
mass and age of the universe are a function of the an- 
gle (f) and the present epoch value of the development 
angle tjh- 

Mg = 2.3 X IO^^Mq ((1 + cost^h) cos^ (/)sin0)~^ 



Buote 



& Tsai 1995), and typically there are merger clues in 



the galaxy and gas distributions. 

In the second scenario, a third object of compara- 
ble mass to the two merging clusters is close enough 
to invalidate the two body analysis. In this case, the 
observed distribution of nearby galaxy clusters should 
be enough to determine whether or not the radial in- 
fall model is appropriate. 

4.2. Merger Mass in A2626 



a;. 



Avlos 
lOOOkm/sy VlMpc 



(4-2) 



, no in8 1' (riH - sin tjh) sin riH , , 

tn = 9.8 X 10 yrs - — ^-^5 cot 

(1 — cos r]H) 

'1000 km/s 

Avlos 



VI Mpc; I, 



Using an approximate age of the universe to deter- 
mine rjH, we arrive at the total cluster mass function 
Mg ((/)). Fortunately the dependence of AIg on tn is 
rather weak (see Fig. If the line of sight sepa- 
ration between the subclusters AIlos is measured, (f) 
is determined, yielding the total system mass (which 



To calculate Mg we require (see Eqn. 4-2) the 
observed rest frame velocity difference and the pro- 
jected separation. The rest frame velocity difference 
is Av = 2,486± 112 km/s. The projected separation 
between the centers of mass is assumed to be the dis- 
tance between the bright, central elliptical in group A 
(its position is consistent with the peak in the X-ray 
emission) and the centroid of the 30 group B members 
with measured redshifts (using the centroid of group 
A members slightly increases the projected separa- 
tion and the merger mass). At a cluster distance of 
I = 175h^^ Mpc, the projected separation between 
the two groups is Al± = (0.707 ± 0.052)ft.-i Mpc. 
The uncertainty reflects the statistical uncertainties 
associated with centroiding group B. The fractional 
uncertainty in the merger mass Mg due to Avlos 
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and Al± is 12%. 

Fig. H displays the ratio of the cluster merger mass 
Mq to the sum of the virial masses of groups A and 
B {Ma + Mb = 9.1 x IO^'^H-'^Mq; see Table |l| and 
MGW96) versus the line of sight separation between 
the subclusters, AIlos = Al±/tan(j). The dotted 
line corresponds to tn = 18 Gyr and the solid line 
corresponds to tn = 13 Gyr. The FP distances to 
galaxies within groups A and B yield a ratio of the 
group B distance to the group A distance [Ib/Ia) 
which can be used to calculate Al^os- Specifically, 
we estimate the line of sight separation between the 
subclusters as Al^os ^ ^75h^^{l — Ib/Ia) Mpc. The 
sample of 24 FP distances to groups A and B yields 
A/los = (14± 18)/i-i Mpc. 

Because distance uncertainties in this case are 
large, our method does not yield tight mass con- 
straints in Abell 2626; however, if groups A and B are 
merging (99% confidence from FP distances and 93% 
confidence from galaxy magnitude distributions), it is 
clear that the virial masses of groups A and B under- 
estimate the total system binding mass. Specifically, 
the sum of the masses of groups A and B must be 
at least 1.5 x lO^^M©, compared to the virial sum of 
9.1^2 7 X lO^Mo (MGW96), a factor of 1.65 higher. 
The virial mass range corresponds to 90% statistical 
confidence limits, where it is assumed the mass uncer- 
tainties are dominated by the velocity dispersion un- 



certainties (Heisler, Tremaine & Bahcall 1985). The 
virial sum is inconsistent with the minimum merger 
mass at the ~ 3a level. 

5. Discussion 

The study of double clusters with the fundamental 
plane holds promise for significant progress in two ar- 
eas. First, the study of appropriately isolated nearby 
double clusters should lead to new constraints on the 
total cluster binding mass. These constraints will not 
only provide an independent test of the virial, hy- 
drostatic, and weak lensing mass estimators, but will 
also provide information about the mass (and typical 
M/L) outside the virialized cluster region. Second, 
the study of more distant double clusters (at z > 0.1) 
provides a means of directly probing for environmen- 
tally driven biases in FP distance estimates. 

5.1. Fundamental Plane Analysis 

Here we apply this approach to the double cluster 
A2626 at cz ~ 17,500 km/s. Because of the cluster 



distance, it is critical to make seeing corrections while 
extracting 9e and (fJ-e)- Our method is similar to the 
one described by Sagha et al. (1993b), but includes 
a seeing-correction uncertainty which makes the pho- 
tometric parameters less sensitive to the in/exclusion 
of central points in the profile. We also note that 
the cosmological dimming correction to the surface 
brightness has a peculiar velocity dependence which 
may be important in cases where many FP distance 
estimates are combined to produce a single, more ac- 
curate cluster distance (§2.2). 

Through numerous cross checks we demonstrate 
that our estimates of the combination log (6*6) — j3 (/ig) 
which appears in the FP are sufhciently accurate that 
their errors make no significant contribution to the FP 
scatter; using a Monte Carlo approach we estimate 
the velocity dispersion uncertainties (§2.3). These pa- 
rameters for the 8 galaxies in group B and the 16 in 
group A are listed in Table ||. 

We combine galaxies from both clusters in a single 
fit, weighting to account for incompleteness (§3.1). 
Our best fit FP coefficients (a = 1.30 ± 0.36 and 
/9 = 0.31±0.06) are statistically consistent with mea- 
surements in other clusters (see Fig. ||). The RMS 
scatter about this plane in Alogi?e is 0.09, corre- 
sponding to 23% distance uncertainties per galaxy. 

We examine the distribution of Mgi, equivalent 
widths with the galaxy spectra from both subclus- 
ters (§3.2). We find that there is no evidence for sig- 
nificant differences in the Mg^-a relation (see Fig. 
^). Our result does not contradict the recently 
noted correlation between cluster environment and 



linestrength (JFK96, |j0rgensen 1997| ); the expected 
offset in Mgi, (given the velocity dispersion differences 
between groups A and B) is too small to detect with 
our data. 

The difference in the fundamental plane zeropoints 
for the two clusters is 7a — 7b = — 0.037 ± 0.046, con- 
sistent with no offset (§3.3). Under the assumption 
that this zeropoint difference is indicative of distance 
differences, log(L»B/L»A) ^ -0.037±0.046, where the 
uncertainty follows from scatter around the median 
distance in each cluster with the best fit FP coeffi- 
cients. We use the 10"* bootstrap resampling simula- 
tions to measure the variation of the distance offset 
as FP coefficients vary. The relative distance con- 
straint is robust enough to reject the Hubble flow hy- 
pothesis {\og{DBlDA) = 0.065) in 99% of the simu- 
lations. Under the Hubble flow hypothesis, the two 
clusters are not interacting gravitationally, and the 
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velocity difFerence reflects pure Hubble flow. An anal- 
ysis of the R band magnitude distributions of con- 
firmed members of both clusters rules out Hubble flow 
with 93% confidence (MGW96). 

If the subclusters were on an outgoing trajectory 
(bound or unbound), then the line of sight distance 
difference would have to be greater than their cur- 
ren t recession velocitv times the age of the universe 
r n? > vr.n.^tH ^ 34 Mpc) which is greater than 



the Hubble fiow value: \og{DB/DA) > 0.0645; thus, 
models where groups A and B are on outgoing tra- 
jectories are ruled out with higher confidence than 
pure Hubble fiow. Therefore, our relative distances 
indicate that the subclusters are bound and infalling 
with 99% confidence (see Fig. ^). 

5.2. Radial Infall Model 

Assuming bound and infalling subclusters, we model 



these kinematics with a radial infall model (Beers 



Seller fc Huchra 1982 ); this two body model is appro- 
priate for separations large compared to the scale of 
the virialized region of the cluster and in cases where 
there are no other massive clusters in the neighbor- 
hood. Given the projected separation, line of sight 
velocity difference, and line of sight separation, the 
radial infall model provides an estimate of the total 
cluster binding mass which depends weakly on the 
adopted age of the universe. From Fig. ^, the total 
binding mass is at least 1.65 times higher than the 
virial sum, a difference significant at ~ 3cr. The total 
binding mass could be much larger than this mini- 
mum. 

The differences in these two mass estimates could 
indicate (1) deviations from radial infall, (2) an error 
in the virial estimate, or (3) significant mass beyond 
the virialized region. The X-ray image of A2626 pro- 
vides no indication of interactions between the two 
subclusters, and there is no third subcluster, so the 
radial infall model should be appropriate in A2626. 
If the virial estimators are biased low, then the radial 
infall model indicates that AI/Ln ~ lOOO/i (see Ta- 
ble which would be a surprisingly high value (e.g. 
Ramella et al. 1989). These arguments suggest that 
the minimum merger mass is larger than the virial 
mass because it is sensitive to mass in the cluster in- 
fall region. 

Taking the merger mass as a lower limit on the 
mass within the cluster infall region, we can address 
whether there is evidence for variation of the mass to 



light ratio between the cluster core and infall region. 
The galaxy light is calculated within a region which 
only extends to a projected radius of 1.5/i^^ Mpc; 
presumably both the galaxy and mass distributions 
extend beyond this region. We have no CCD pho- 
tometry over the larger region to directly measure 
the total cluster light. Assuming that the projected 
galaxy distribution falls off as ~ R~^-^ (e.g. Mohr 
et al. 1996| ) with a similar luminosity function to 



the galaxies in the central region, the total cluster 
light grows as ~ ^fR for R larger than the core ra- 
dius. Under these assumptions, within a diameter of 
^ 4/i^^ Mpc there is enough light to yield a mass-to- 
light ratio of M/Lr ~ 600/i throughout the cluster 
even with the higher mass estimate from the radial 
infall model. Of course, the binding mass could be 
many times larger than the minimum value, so the 
mass-to-light ratio could increase significantly out- 
side the virialized cluster core. A study of the cluster 
light extending to larger radius and a larger sample of 
FP distances would yield more concrete information 
regarding any variation in the cluster mass-to-light 
ratio in A2626. 

FP studies of other double clusters, coupled with 
(1) a more detailed accounting of the radial distribu- 
tion of cluster light and (2) larger redshift samples to 
enable a more accurate estimate of the virial mass, 
should provide interesting constraints on the distri- 
bution of cluster mass and possible differences in the 
efficiency of galaxy formation beyond the virialized 
region. 

We thank Roberto Saglia and Glenn Baggley for 
their comments on an earlier version of this manuscript 
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our measurements to their unpublished results. We 
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attempting additional imaging. This research was 
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A. Effects of Culling the Galaxy Sample 



Four of our 24 galaxies show peculiar morphologies; here we detail the peculiarities and demonstrate that 
eliminating them from the analysis does not significantly alter our conclusions. The following are morphological 
types which we derived by comparing deep 2.4 m images and the Carnegie Atlas of Galaxies ( ^andage fc Bedke 
1991): 



• AH: bright nucleus with possible bar; probable SBO 

• AL: possible spiral structure; SBO or Sa 

• AM: bright nucleus with possible bar; probable SBO 

• AP: appears to have a dust ring; probable RSBO 



Without these four galaxies, Group A (B) contains 12 (8) galaxies; we now apply the analysis detailed in §3 to 
this smaller sample. The best fit form of the FP is 

log 6ie = (1.39 ± 0.65) logo- + (0.32 ± 0.06) (/ie) + 7, (Al) 

and the scatter around this plane in logi?e is 0.080 (20% distance uncertainty). These coefficients are consistent 
with the values from the full sample, but the scatter is somewhat smaller and the a uncertainty is larger. The FP 
coefficient uncertainties are determined (as before) through bootstrap resampling and refitting the sample. 

The RMS scatter around the median distance estimator in group A (B) is 0.066 (0.115); the zeropoint difference 
corrected for homogeneous Malmquist bias is 

A7 = -0.026 ±0.045. (A2) 
The zeropoint difference and uncertainty (from the scatter) are somewhat smaller than for the whole sample (see 



Eqn. 3-4). This measurement is inconsistent with the Hubble flow hypothesis at ~ 2. Oct or 97.7% confidence. The 
bootstrap refitting produces a distribution of A7 which includes the variation in the FP coefficients; 5000 refits of 
this smaller sample rule out the Hubble flow hypothesis with 98.2% confidence. In summary, removing 4 galaxies 
which have morphological peculiarities slightly reduces the signiflcance with which the Hubble flow model can be 
rejected in Abell 2626. 



Table 1 
subcluster properties 



ID 


N 


V 


CT 


M„,[1014Mo] 




A 
B 


67 
30 


16,533±141 
19,164±138 


658^^}^ 


9 q+1.4 
^■'-'-0.7 


630 
570 



Intervals are statistical 90% confidence limits 
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Table 2 
Galaxy Properties 



Tag 


RA(1950) 


Decl 




cz 


(7 


Co- 


Mr 


(Me) 


( 


9° 

e 






AA 


23 


34 


29, 


.30 


20 


22 


22. 


.3 


17,538 


287 


12 


14, 


.06 


19, 


.79 


6, 


.232 


4. 


,83 


0.53 


AB 


zo 


oo 


1 8 

lo. 


91 


90 




ou. 


.7 


17,159 


223 


22 


lO, 




1 Q 




5, 


.477 


A 


,ou 


0.34 


AC 


23 


34 


22, 


.71 


20 


56 


55. 


.6 


16,819 


218 


22 


14, 


.37 


19, 


.56 


4, 


.846 


4. 


,41 


0.79 


AD 


23 


33 


40, 


.11 


20 


45 


31, 


.4 


17,911 


190 


oo 

11 


14, 


.75 


19, 


.93 


4, 


.817 


4. 


,64 


0.63 


AE 


23 


33 


23, 


.09 


20 


42 


02, 


.1 


15,941 


205 


28 


15, 


.03 


19, 


.22 


3 


.063 


4, 


,42 


0.55 


AF 


23 


35 


33, 


.91 


21 


07 


54, 


.7 


16,728 


191 


12 


14, 


.97 


19, 


.23 


3, 


.161 


4. 


,17 


0.74 


AG 


23 


32 


40, 


.07 


21 


00 


01, 


.6 


15,979 


163 


19 


14, 


.71 


19, 


.33 


3 


.727 


3. 


,83 


0.85 


AH 


23 


33 


40, 


.68 


20 


42 


35, 


.6 


16,644 


161 


6 


15, 


.02 


20, 


.15 


4, 


.732 


4, 


,54 


0.50 


AI 


23 


34 


32, 


.46 


20 


18 


14, 


.5 


17,804 


185 


7 


14, 


.93 


19, 


.24 


3, 


.231 


4. 


,14 


0.84 


AJ 


23 


34 


35, 


.93 


20 


59 


41, 


.0 


16,665 


178 


35 


15, 


.35 


18, 


.94 


2, 


.319 


4. 


,67 


0.90 


AK 


23 


33 


58, 


.71 


21 


01 


12, 


.4 


16,856 


202 


20 


15, 


.35 


19, 


.23 


2, 


.564 


4, 


,86 


0.70 


AL 


23 


34 


08, 


.75 


20 


49 


29, 


.9 


17,645 


152 


3 


14, 


.61 


19, 


.67 


4, 


.554 


3. 


,67 


0.39 


AM 


23 


33 


49, 


.52 


20 


49 


12, 


.1 


16,129 


252 


21 


14, 


.40 


21, 


.03 


9, 


.403 


4. 


,87 


0.64 


AO 


23 


33 


57, 


.86 


20 


51 


11, 


.8 


16,113 


190 


27 


15, 


.91 


17, 


.66 


0, 


.997 


4, 


,68 


1.00 


AP 


23 


33 


47, 


.81 


20 


49 


21, 


.8 


16,859 


172 


8 


15, 


.07 


19, 


.70 


3, 


.743 


4. 


,86 


0.30 


AZ 


23 


35 


19, 


.60 


20 


30 


15, 


.3 


16,063 


118 


11 


16, 


.02 


19, 


.74 


2, 


.136 


3. 


,89 


0.70 


BA 


23 


34 


07, 


.96 


20 


31 


29, 


.6 


19,729 


148 


3 


14, 


.84 


19, 


.19 


3 


.351 


3, 


,80 


0.69 


BB 


23 


33 


19, 


.80 


20 


33 


12, 


.0 


19,369 


219 


21 


14, 


.31 


19, 


.61 


4, 


.900 


4. 


,27 


0.74 


BC 


23 


33 


57, 


.29 


20 


35 


11, 


.2 


19,240 


187 


17 


14, 


.36 


20, 


.06 


6, 


.243 


4. 


,42 


0.69 


BD 


23 


33 


08, 


.31 


20 


48 


02, 


.8 


19,150 


214 


10 


14, 


.21 


20, 


.48 


8, 


.148 


4, 


,45 


1.00 


BE 


23 


33 


17, 


.67 


20 


22 


54, 


.2 


19,215 


133 


3 


15, 


.18 


19, 


.17 


2, 


.836 


4. 


,27 


0.62 


BF 


23 


33 


33, 


.61 


20 


49 


36, 


.5 


19,825 


119 


3 


15, 


.89 


19, 


.77 


2, 


.705 


3. 


,79 


0.38 


BG 


23 


33 


17, 


.50 


20 


52 


03, 


.1 


18,818 


133 


6 


16, 


.02 


19, 


.71 


2, 


.431 


4, 


,61 


0.68 


BH 


23 


33 


12, 


.20 


20 


35 


20, 


.1 


19,136 


151 


5 


17, 


.09 


18, 


.34 


0, 


.761 


4. 


,98 


0.33 



'^in arc seconds 

^bulge luminosity fraction 
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1.4x10* 1.6x10* 1.8x10* 2x10* 2.2x10* 
v[km/ s] 



Fig. 1. — The galaxy velocity distribution around Abell 2626. The histogram contains 108 galaxy velocities from the 
redshift survey of MGW96. Group A is centered at t; = 16, 500 km/s, and group B is centered at u = 19, 150 km/s. 
The high fraction of galaxies with emission lines in the group of 11 centered at ?; ^ 21,500 km/s (73% compared 
to 39%/33% for groups A/B) indicates they are part of a low density background structure. 
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Fig. 2. — Two projections of the FP. Both plots contain group A (sohd) and B (hollow) galaxies. The lines define 
the best fit group A (solid) and B (dashed) FP; offsets in the fits reflect subcluster distance differences. The best fit 
FP (with weighting to account for incompleteness- see text) is of the form log 9e = 1.30 log a + 0.31 (/Xe) + 7cj where 
7c; is the distance dependent cluster zeropoint. The RMS scatter around the best fit in log^e is 0.09, corresponding 
to a 23% distance uncertainty per galaxy. 
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Fig. 3. — Plot of log (Mgi,) vs log (a) for the galaxies in groups A (solid) and B (hollow). The lines mark the best fit 
relations for groups A (solid) and B (dashed) constrained to have the slope found using a much larger galaxy sample 
(J0rgensen 1997). The apparent zeropoint offset between the two samples is statistically insignificant (0.012 ±0.019; 
see text), providing no evidence for significant stellar population differences between the subclusters. 



17 



2x10* 



>, 1.8x10" 
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<Vb> o 



<Va> 



3 

SG 
O 



-0.2 



-0.1 



0.1 



0.2 



log (D/DJ 



Fig. 4. — Plot of fundamental plane relative distances versus redshift for galaxies in groups A (solid) and B 
(hollow). Galaxies in groups A and B are clearly offset in velocity (see Fig. |l|), and also appear to be offset in 
distance. Vertical lines mark the (homogeneous Malmquist bias corrected) median distances to groups A (solid) and 
B (dashed). Horizontal lines mark the mean velocities of groups A (solid) and B (dashed). The large star marked 
with Dnubbio lies at the expected position of group B if its trajectory relative to group A were pure Hubble flow. 
The median distances to groups A and B indicate that log (Db/Da) = — 0.037 ± 0.046; the Icr distance uncertainty 
is indicated by the error bar on the vertical, dashed line. The sample of 24 galaxies rules out the Hubble flow 
hypothesis, Db = Dnubbio and log (Db/Da) — 0.065, with 99% confidence (see §3.3). 
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Fig. 5. — The ratio of the cluster merger mass Mq to the sum of the group A and B virial masses Ma + Mb = 
9.1 X lQ^'^h~^ Mq as a function of the hue of sight separation AIlqs of the two groups. The two solutions correspond 
to different assumptions regarding the age of the universe tn- tn = 13 Gyr (solid) and tn = 18 Gyr (dotted). The 
statistical uncertainty in the merger mass Mq is 12%. If groups A and B are merging, as suggested by the FP 
distances and the galaxy magnitude distributions, the merger mass is at least 1.65 times greater than the sum of 
the virial masses. The FP distances indicate that AIlqs = 14 ± 18 Mpc. 
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